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We report the preliminary results of the studies of hadron spectrum under the background of abelian and 
monopole gauge fields in quenched Wilson SU(3) QCD. Abelian gauge fields alone reproduce the same chiral limit 
as in the full case. Critical hopping parameter k c and m p are the same in both cases. We need more time to get 
a definite result in the case of monopole background. The photon contribution do not produce any mass gap in 
the chiral limit (k = k c ~ 0.17). The behavior is similar to those in the free photon case for k c = 0.125. 



1. Introduction 

Recent Monte-Carlo simulations suggest that 
abelian monopoles after abelian projection are a 
key quantity of confinement mechanism in QCD 
[1-3]. Especially the abelian projection in the 
maxumally abelian (MA) gauge is interesting. 
Then mass generation from scale invariant QCD 
may be explained by the abelian monopoles alone. 
The expectation is supported by the preliminary 
data of recent simulations in quenched SU (2) 
QCD with Kogut-Susskind fermions[4]. The pur- 
pose of this study is to report preliminary results 
of the studies of hadron spectrum under the back- 
ground of abelian and monopole gauge fields in 
quenched SU(3) QCD with Wilson fermions. 

2. Maximally abelian gauge in SU(3) QCD 

We adopt the usual SU(3) Wilson action. The 
MA gauge is given by performing a gauge trans- 
formation, U(s,fi) = V(s)U(s,ii)V~ 1 (s+jj), such 
that 

R = E (\Un{ 8 ^)? + \U22(s,fi)\ 2 + \U 33 (s,fi)\ 2 ) 

S ,/J, 

is maximized. Then a matrix 

+U\s - fi,n)k a U{s - fi,n),k a ] 
is diagonalized.Here 

Ai = diag(l, -1, 0), A 2 = diag(-l, 0, 1), 



A 3 = diag(0,l,-1). 

After the gauge fixing is over, we can extract an 
abelian link field, 

U(s,n) = A(s,n)u(s,n), 

where 

u( S ,fi) = dmg(e* s(1) ( s ^,e* s(2) ( s ^,e* s(3) ( s ^) 
<i>( s ,V) = X! ar S(^»( s '^)) lmod2^ • 

i 

u(s, /i) is a diagonal abelian gauge field. 

3. Abelian gauge fields from monopole 
Dirac strings 

A plaquette variable is given by fjjj(s) = 
d^\s)-d v e ( p(s), where 

s' 

and D(s — s') is the lattice Coulomb propagator. 
Define abelian gauge filed as 

e ( p( s ) = -Y J D ( s -^) d 'J^)- 

s' 

The gauge fields automatically satisfy the Landau 
gauge condition. 

Extract the Dirac strings from the field 
strength, satisfying nfy(s) = 0: 

/#(*) = /$W + 2™«(«), 
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Then an abelian gauge field from Dirac string is 
S^BW = -2^D(«- 8 ')fl:n«( S '), 

s' 

whereas an abelian gauge field from the photon 
is 

C ( ^) = -E^( s - s W«(s')- 

s' 

Abelian link fields from Dirac string and pho- 
ton are constructed as 

u Ds ( S ,v) = diag(e< ,(1) ('),e< ,(a) ('),e< ,(3) (')), 
u p \ 8 ,ri = diag(e i9 r (1) ('),e i9 r (a) ('),e i9 r (3) (')). 

4. Inverse of Wilson quark matrix 

Now let us evaluate the inverse of Wilson 
fermion matrix 

+ (1 + JnWH-s - fi,fi)6 s -ji !S >] 

for 

• U{s,n) (full) 

• u(s,/i) (abelian) 

• u Ds (s,/j) (Dirac string) 

• u ph (s,/j) (photon), 

and then measure hadron mass spectrum. 
The simulations are performed as follows: 

• Our calculation is usually done on a 8PE's 
system of a vector parallel supercomputer 
Fujitsu VPP500. Lattice size adopted is 
16 3 x 32. /3= 5.7. 

• Antiperiodic (Periodic) boundary condi- 
tions in the time (space) direction. 

• Number of sweeps for thermalization is 
3000. 

• Number of sweeps to get independent con- 
figurations is 1000. 

• Number of configurations for average is 20. 



• The criterion of the MA gauge condition 
is I °ff diagonal part of X\ < 10 -7 where 
(X is the operator to be diagonalized.) 
Over-relaxation method with u> = 1.93 is 
used. 

• The criterion of the solution of Wilson ma- 
trix inversion is |residue| < 10 -5 . The 
hopping parameter expansion and the CG 
method with red-black preconditioning is 
used. 

• The mass fitting function is ccosh(m(t — 
Nf/2)) and the mass fitting range is 4 < 
N t < 28. 

The CPU time , the number of terms of the 
hopping parameter expansion and the number of 
CG iterations in a quenched propagator compu- 
tation are listed in Table 1. 



5(7(3): 


K 


time(sec) 


HPE 


CG 


0.160 


230 


1200 


240 


0.165 


551(16PE) 


1200 


3000 


abelian: 


K 


time(sec) 


HPE 


CG 


0.150 


482 


360 


1560 


0.155 


910 


300 


3600 


0.1575 


1188 


120 


4440 


0.160 


2074 


300 


8040 


0.1625 


3000(16PE) 


300 


20400 


monopole: 


K 


time(sec) 


HPE 


CG 


0.135 


345 


1200 


600 


0.140 


429 


180 


1560 


0.1425 


346(16PE) 


120 


2160 



Table 1 

Various parameters of the simulations. 



5. Results 

Examples of hadron propagators with abelian 
and monopole backgrounds are shown in Fig. 1 
and Fig. 2. 
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abelian propagator ( k = 0.150) 
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Figure 1. Hadron propagators with abelian back- 
ground for k = 0.150. 

The estimated masses of w, p and nucleon are 
listed in Table 2 [5]. 

The squared masses of 7r and p and the mass 
of nucleon and p are plotted versus 1/k in Fig. 3 
and Fig. 4. Abelian gauge fields alone reproduce 
the same chiral limit as in the full case. Critical 
hopping parameter k c and m p are the same in 
both cases as expected. The abelian case gives a 
smaller ratio m w /m p at the same k. At present we 
have obtained the ratio as small as m^/m p ~ 0.48 
at k = 0.1625. 

The abelian case gives a smaller ratio m p /m p 
at the same k. Namely the 0(a) correction of this 
case is smaller in the abelian case. At present we 
have obtained the ratio as small as m p /m p ~ 1.32 
for m p a = 0.55 ( at k — 0.1575). 

We need more time to get the monopole con- 
tribution clearly and it is in progress. 

Fig. 5 shows the squared masses of 7r and p ver- 
sus 1/k under the photon and free backgrounds. 
The photon contribution do not produce any 
mass gap in the chiral limit (k = k c ~ 0.17). 
The behavior is similar to those in the free pho- 
ton case, although k c = 0.125 in the free case. 

In summary, the above preliminary data sug- 
gest abelian and also monopole gauge fields are 
essential for the mass generation of hadrons in 
QCD. More data are needed to get a definite con- 
clusion and it is under progress. However the 



monopole propagator ( k = 0.140) 
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Figure 2. Hadron propagators with monopole 
background for k = 0.140. 



data are consistent with the results of recent sim- 
ulations in quenched SU (2) QCD with Kogut- 
Susskind fermions[4]. 
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5(7(3): 


K 


m w 


m p 


m N 


0.160 
0.165 
0.169(0 


0.689(3) 
0.459(4) 
0.0 


0.813(4) 
0.676(1) 
0.543 


1.256(18) 
1.045(25) 


abelian: 


K 


m w 


m p 


m N 


0.150 

0.155 

0.1575 

0.160 

0.1625 

0.170 (O 


0.545(6) 
0.442(5) 
0.364(4) 
0.300(4) 
0.244(4) 
0.0 


0.630(5) 
0.543(4) 
0.514(6) 
0.492(9) 
0.503(23) 
0.553 


1.00(2) 
0.81(3) 
0.68(7) 


monopole: 


K 


m w 


m p 


m N 


0.135 

0.140 

0.1425 

0.150 

0.1575 


0.671(5) 
0.533(3) 
0.500(5) 
0.314(6) 
0.139(18) 


0.713(5) 
0.606(3) 
0.581(6) 
0.442(14) 
0.247(102) 


1.29(7) 
1.04(3) 
0.91(5) 



Table 2 

Hadron masses. 
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Figure 3. The squared masses of 7r and p versus 
1/k. 
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Figure 4. The masses of nucleon and p versus 
1/k. 
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Figure 5. The masses of nucleon and p versus 
1/k. 
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